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Bose-Einstein correlations of charged kaons are used to probe Au+Au collisions at ys NN = 200 
GeV and are compared to charged pion probes, which have a larger hadronic scattering cross section. 
Three dimensional Gaussian source radii are extracted, along with a one-dimensional kaon emission 
source function. The centrality dependences of the three Gaussian radii are well described by a single 
linear function of N^ t with zero intercept. Imaging analysis shows a deviation from a Gaussian 
tail at r > 10 fm, although the bulk emission at lower radius is well-described by a Gaussian. The 
presence of a non-Gaussian tail in the kaon source reaffirms that the particle emission region in a 
heavy ion collision is extended, and that similar measurements with pions are not solely due to the 
decay of long-lived resonances. 
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Experiments at the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory have re- 
vealed that collisions of Au ions at y/s^ = 200 GeV 
produce a new form of matter which is opaque to jets and 
exhibits anisotropic flow consistent with perfect fluid hy- 
drodynamics [3, 3] • Studies of the space-time evolution 
of the collisions are needed to elucidate the properties of 
the hot, dense, and strongly interacting matter, probe the 
time scale and degree of thermalization, and investigate 
the order of the deconfinement phase transition. 

Two-particle interferometry, also known as HBT af- 
ter the radio astronomers R. Hanbury Brown and R.Q. 
Twiss (3, is a powerful tool for measuring the space- 
time extent of particle-emitting sources. In elementary 
particle and nuclear physics, enhanced production of like- 
sign pions with small relative momenta was discovered 
experimentally and explained by the Bose-Einstein sym- 
metrization of identical bosons [4| . Correlations are pro- 
duced by the combination of quantum mechanical inter- 
ference of identical particles and strong and/or electro- 
magnetic final state interactions such as Coulomb re- 
pulsion for same-sign charged pairs. HBT radii refer 
to Gaussian measures of source sizes on the femtometer 
scale. 

Although the traditional HBT analyses are constrained 
by the assumption of a Gaussian distribution of particle 
emission, recent detailed measurements of pion emission 
sources using an imaging technique show a non-Gaussian 
structure for the two-particle source region above ~20 
fm Q , suggesting the possibility that decays of long-lived 
resonances or a temporal component of the source con- 
tribute to the non-Gaussian tail [f|. While charged pi- 
ons are strongly affected by rescattering among hadrons 
and decays of hadronic resonances, charged kaons have 
smaller rescattering cross sections than charged pions 
{&k-n < °V-tv) and are less affected by resonance de- 
cays. Until recently, no full hydrodynamic calculation 
has accurately predicted particle spectra and HBT radii 
for pions, and none can simultaneously describe the mo- 
mentum asymmetry measurements of flow. The measure- 
ment of the kaon source describe herein add an important 
new constraint to address this "HBT puzzle" @, S|- 

An angle-averaged one-dimensional Gaussian measure- 
ment of correlations of neutral kaons by STAR 0] sug- 
gests that the transverse mass dependence for neutral 
kaons and charged pions falls on one universal curve. In 
this paper, 3D Gaussian HBT correlations of like-sign 
kaons are presented in three transverse momentum bins 
for 0.3 < pt < 1.5 GeV/c and three collision centrality 
bins. The resulting HBT radius parameters for kaons are 
compared to those of like-sign pion pairs L 10J. In addi- 
tion we present ID emission source functions for charged 
kaons in relativistic heavy-ion collisions. 



This analysis of 2004 data from the PHENIX detec- 
tor [ll[ uses ~ 600 million minimum bias events, which 
are triggered by the coincidence of the Beam-Beam Coun- 
ters (BBC) and Zero-Degree Calorimeters (ZDC) with 
collision vertex \z\ < 30 cm. A Monte Carlo Glauber 
model [HI is used to match the observed BBC and ZDC 
distributions and to bin the data according to the number 
of nucleons participating in the collisions, -/Vp art . 

Charged kaons are tracked and identified using the 
drift chamber (DC), pad chambers (PC1,PC3) and PbSc 
Electromagnetic Calorimeters (EMCal) to cover pseu- 
dorapidity \r/\ < 0.35 and azimuthal angle A<f> = n/2 
(A<f> — 7r/4) in one (and the other) central arm. A 
track model provides a 3-dimensional trajectory and mo- 
mentum vector for charged particles based on DC and 
PCI information with a momentum resolution of Sp/p ~ 
0.7% © 1.0% x p (GeV/c). Backgrounds are reduced by 
requiring 2 a position match between track projections 
and EMCal hits, and 3 a match for PC3. Kaons are 
separated from pions up to pr ~0.9 GeV/c using timing 
information from BBC and EMC. Particles at higher px 
that fall within 2 a of the ideal mass-squared for kaons 
but > 3 a away from the peak for pions or (anti-)protons 
are identified as kaons. The contamination level is ~4% 
from pions, and ~1% from protons at pr ~ 1-5 GeV/c. 

The correlation function is experimentally measured 
as C2(q) = A(q)/B(q) where A(q) is the relative mo- 
mentum (q) distribution of actual pairs obtained by all 
possible combinations of pairs within the same events 
and B(q) is the background pair distribution from mixed 
events. Two-track detection inefficiencies for charged 
kaons that traverse the DC and EMCal in close proximity 
have been carefully studied with Monte-Carlo detector 
simulation and the actual pair distribution is corrected 
by the MC efficiency factors. After pair selection cuts to 
remove track splitting and merging (see [lfj| for details), 
~15 million positive kaon pairs and 14 million negative 
kaon pairs remain. 

To measure multi-dimension source sizes, q is decom- 
posed into standard "side-out-long" axes [13j: for which 
<7iong is parallel to the beam axis, q ou t is parallel to the 
transverse momentum of the pair (fcr = (pit +P2t)/2), 
and g s ide is orthogonal to both gWg and g ut • This anal- 
ysis is performed in the Longitudinally Co-Moving Sys- 
tem (LCMS) defined as p±z — ~Ptl- For the treatment 
of charged kaons emitted away from the central region 
(core), we adopt an effective core-halo Coulomb correc- 
tion, proposed by Bowler and Sinyukov 14| , in which the 
3D Gaussian fit function is given by 



Qcore 



Cr° = [A(l + G)] F C + [1 - A] , (1) 



where the Coulomb correlation function Fc is iteratively 
evaluated from the Coulomb wave function of kaon pairs 
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FIG. 1: (color online). 3D correlation function of charged 
kaon pairs measured for 0.3 < &t < 1.5 GeV/c at — 30% 
centrality in Au+Au collisions at ^/s NN = 200 GeV. The curve 
is a 3D function is a Coulomb uncorrected fit with Eq. [T]and 
sliced over the lowest 40 MeV/c in the orthogonal directions. 



assuming a spherical Gaussian source. The Gaussian cor- 
relation function in the side-out-long decomposition is 
determined by 



G = cxp [-R 2 sidc q : 



Jsidc 



d2 2 

^outyout 
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out 



(2) 



The systematic error estimate incorporates a contri- 
bution from the Coulomb interaction of the source halo 
using a prescription developed by Maj and Mrowczyn- 
ski [l5(. The fitted i? S idc and i?i on g are Gaussian mea- 
sures of the spatial lengths of homogeneity, where parti- 
cles of similar momenta are emitted [Iff , in the transverse 
and longitudinal directions at freeze-out. i? ut contains 
a contribution from the duration of the particle emission 
in addition to the spatial length [l7T |. Note that an out- 
long cross-term vanishes in the expression for G for our 
\r]\ < 0.35 acceptance at midrapidity [HI]. The fitted A is 
empirically defined and includes contributions from mis- 
identified particles (1 — f) 2 along with components of the 
source that are not well resolved by the Gaussian fit. 

Figure Q] shows the 3D correlation function of charged 
kaons without the Coulomb correction measured for 
0.3 < fc T < 1.5 GeV/c at - 30% centrality in Au+Au 
collisions at ^/s NN = 200 GeV and the fit with Eq. JTJ). 
Separate fits to the 2K + and 2K~ correlation functions 
were performed, yielding consistent results for all &t and 
centrality bins. 

Panels (a)-(e) in Fig. [2] show the HBT radius param- 
eters of charged kaons for 0.3 < fcx < 1.5 GeV/c as 

1/3 

functions of N vaxt , which is proportional to the trans- 
verse radius of the initial collision volume. Similar to 
pions, kaon radii are well described by linear functions of 

1/3 



7V p ^ rt . Because initial fits yielded slopes that were consis- 
tent for all radii and intercepts that were consistent with 
zero, the all three radii were fit to a single linear function 

1/3 

with zero intercept: Ri =pi*Af p ^ rt , withpi = 0.51±0.01 
and x 2 / n df = 14.6/8. We note that similar fits to the 
pion radii yield non-zero intercepts. Although pions and 
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FIG. 2: (color online). 3D Gaussian HBT radius parameters 
extracted from Eq.{T} for charged kaon pairs as a function of 
(left) iVp4 3 t measured for 0.3 < k T < 1.5 GeV/c and (right) 
mi measured for — 30% centrality in Au+Au collisions at 
^/s NN = 200 GeV. In panels (b)-(d) the dashed line is a fit 

with pi * jVp^ 3 t . 



kaons are measured in a similar range (0.2 < fcx < 2.0 
GeV/c for pions), the higher transverse mass (tot = 
\/fcT 2 + to 2 ) for kaons ((m T ) -0.89 GeV/c 2 ) than for 
pions ((tot) —0.47 GeV/c 2 ) leads to smaller radii, as ex- 
pected from (mx) scaling 0, ljj 20l[2lT|. 

Panels (f)-(j) in Fig. [2] show the tot dependence of 
the radius parameters for kaons in 3 different tot bins at 
— 30% centrality, compared with pions [To| . HBT radii 
are quite consistent at the same tot, clearly indicating 
that the radii follow tot scaling. The -Rout/^sidc ratio 
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for kaons is ~ 1.0- 1.2 which is consistent with the value 
for pions at low tut- Kaon HBT results from a 2D+1 
hybrid, hydrodynamic + UrQMD calculation (open cir- 
cles) [22J show slightly larger sidewards radii than the 
data, and the outwards and longitudinal components arc 
too large by a factor of 2-3. A more recent 1D+1 hybrid 
calculation [23] for kaons (open squares), which assumes 
flat rapidity distribution and axial symmetry, compares 
more favorably, matching all radii to within systematic 
and statical errors. This calculation incorporates pre- 
equilibrium flow and a lattice-inspired equation of state, 
which are two features lacking in earlier calculations of 
HBT radii. Although promising, these theoretical results 
remain to be verified with full 3D+1 calculations that can 
also reproduce the elliptic flow. 

Recent femtoscopic measurements 0, Q , which use an 
imaging technique Q revealed that the emission source 
function of charged pions has a non-Gaussian tail which 
cannot be resolved with traditional Gaussian fitting tech- 
niques. In the imaging scheme, the correlation function 
is expressed by the Koonin-Pratt equation [H, [2(| 

C 2 (q)-1= / drK(q,r)S(r), (3) 

where the kernel K(q, r) is the relative wave function 
as |<i>q ^ (r) | 2 — 1 that describes the propagation of pairs 
emitted with relative separation r and relative momen- 
tum q in the Pair Center-of-Mass System (PCMS). S(r) 
is the emission source function of pairs. 

The filled squares in Fig. [3] (a) show the ID kaon corre- 
lation as a function of the invariant relative momentum of 
the pair (gi nv = y/—(pi — p2) 2 /2). The ID source func- 
tion S(r) imaged from C^Qinv) is shown by filled circles 
in Fig. [3](b). 

In this analysis, input parameters that govern the 
imaging procedure [3] were selected to minimize the x 2 
between the data and the restored C2(ginv)(x 2 / n df ~ 1), 
shown by open circles in Fig. [3] (a). The solid curve 
shows the traditional Gaussian source function, obtained 
by angle-averaging the 3D HBT radius parameters (A, 
-Rsidc, ^out: R\ong) m the PCMS frame, the same frame 
in which the imaging is performed. 

The imaged S(r) exhibits a non-Gaussian tail at r > 10 
fm. This excess corresponds to the deficit in the <7i nv ^ 
20 MeV/c region of the angle-averaged Gaussian curve 
of Fig. [3] (a) , and is also visible in the 3D Gaussian slices 
in Fig. [TJ The S(r) for pions in the same fcx range shows 
a similar trend. The deviation from a Gaussian in the 
shape of the S(r) indicates that the particle emission re- 
gion is extended, and a similar non-Gaussian tail in the 
pion source is not solely the result of long lived reso- 
nance decays such as the w, although a less prominent 
contribution from the K* is likely. The observation of a 
more substantial non-Gaussian tail for kaons than for pi- 
ons is qualitatively consistent with a hadronic resonance 




q. nv [MeV/c] 




r[fm] 

FIG. 3: (color online), (a) (filled squares) measured C2(ginv). 
(open circles) restored C2((?inv) from imaged S(r), com- 
pared with (solid curve) angle- averaged Gaussian C2(<jw) for 
charged kaon pairs measured for 0.3 < fcx < 0.9 GeV/c at 
- 30% central Au+Au collisions at v /s^7 = 200 GeV. (b) 
(filled circle) Imaged kaon S(r) compared with (solid curve) 
angle-averaged Gaussian S(r). Error bars are statistical only 
and boxes indicate the total systematic errors, (open triangle) 
S(r) for charged pion pairs for the same kx region. For the 
pion S(r), error bars include both statistical and systematic 
errors. 

cascade model with a time dependent density for an ex- 
panding source, in which the larger mean free path for 
kaons leads to an extended emission region [2?], H3| • De- 
tailed measurements with 3D HBT imaging of kaons, or 
ID imaging of more species probing different hadronic 
cross sections will determine contributions from other ki- 
netic effects to S(r). 

In summary, we have measured Bose-Einstein correla- 
tion functions of charged kaon pairs in Au+Au collisions 
at v /s^7 = 200 GeV. The 3D HBT radii i? sidc and i?i ong 
are consistent for pions and kaons at the same N pai t and 
tot- The ID emission source function for kaons extracted 
by imaging shows a non-Gaussian tail at distances greater 
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than 10 fm. This tail represents a direct measurement of 
the ID length of homogeneity of the particle emission 
source and is not due primarily to resonance decays. 
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